Observations of two-dimensional quantum oscillations and ambipolar transport in the 
topological insulator Bi2Se3 achieved by Cd doping 
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We present a defect-engineering strategy to optimize the transport properties of tlie topological 
insulator Bi2Se3 to show a high bulk resistivity and clear quantum oscillations. Starting with a p- 
type Bi2Se3 obtained by combining Cd doping and a Se-rich crystal-growth condition, we were able 
to observe a p-to-n-type conversion upon gradually increasing the Se vacancies by post annealing. 
With the optimal annealing condition where a high level of compensation is achieved, the resistivity 
exceeds 0.5 ficm at 1.8 K and we observed two-dimensional Shubnikov-de Haas oscillations composed 
of multiple frequencies in magnetic fields below 14 T. 

PACS numbers: 73.25.-(-i, 74.62. Dh, 72.20.My, 73. 20. At 



I. INTRODUCTION 



The three-dimensional (3D) topological insulator (TI) 
realizes a novel quantum state of matter where a non- 
trivial Z2 topology of the wavefunction of the bulk va- 
lence band leads to the emergence of a "topological" 
surface state consisting of helically spin-polarized Dirac 
fermions.— "— The peculiar spin texture of the surface state 
holds promise for novel spintronics and fault-tolerant 
topological quantum computing, so there is a rush of re- 
search to address this surface statei^i^ However, most of 
the known TI materials are poorly insulating in the bulk, 
making it difficult to probe the surface state by transport 
experiments. For example, Bi2Se3 is considered to be a 
promising TI material because it has a relatively large 
(^0.3 eV) bulk band gap and a nearly perfect Dirac cone 
as its topological surface statej^i^ however, no matter 
whether it is in the form of bulk crystalj^!^ nanoribbonjS 
or epitaxial thin film,^*^ Bi2Se3 always accompanies a lot 
of Se vacancies (usually ^10^^ cm~^) that act as electron 
donors, and as a result, the residual bulk carriers hinder 
the transport studies of the surface state of this material. 

To achieve a bulk- insulating state in Bi2Se3, doping 
holes to compensate for the residual electrons is a vi- 
able strategy. While this was done through low-level sub- 
stitution of Ca^+ for Bi''+,^i the resulting disorder was 
so strong that no Shubnikov-de Haas (SdH) oscillation 
from the surface state was observed^^ in Bi2_a;Caa;Se3. 
A different strategy was to partially substitute Sb for 
Bi, which apparently reduces the Se vacancies; indeed, 
with a relatively large (~12%) Sb substitution, surface 
SdH oscillations were successfully observed in n-type 
Bi2_a;Sba:Se3, but a very high magnetic field ('^60 T) was 
required for the observationji^ It is to be noted that with 
the Sb doping one can never cross the band gap to reach 
the p-type regime, and hence the tuning of the chemical 
potential to the Dirac point is impossible. This is a pity, 
because Bi2Se3 is attractive for its isolation of the Dirac 
point from the bulk bands. Therefore, it is desirable to 
find a suitable p-type dopant to access the Dirac point 
while keeping the mobility to be sufficiently high for the 



surface state to be studied by the SdH oscillations. 

In this paper, we show that tactful defect engineer- 
ing in Bi2Se3 employing Cd doping in combination with 
Se-vacancy tuning provides a useful means to control 
the chemical potential across the band gap. In the 
literatures^ whereas Cd in Bi2Se3 was shown to be- 
have as an acceptor, Cd-doped Bi2Se3 crystals always 
remained n-type due to the low solubility of Cd atoms 
in Bi2Se3; however, it has been elucidatedi^"— that in- 
creasing the Se content in the Bi-Se melt for the crystal 
growth can suppress the formation of Se vacancies and 
greatly reduce the residual bulk carrier density to the 
level of ^10^^ cm^'^. Therefore, even though the solu- 
bility of Cd is low, one could achieve a p-type behavior 
by combining Cd doping and a Se-rich growth condition. 
Actually, we obtained a p-type sample with this strat- 
egy and, furthermore, starting from the p-type sample, 
we could gradually increase the Se vacancies by careful 
post-annealing and achieve a high level of compensation, 
at which the sample becomes optimally bulk-insulating 
and presents two-dimensional (2D) SdH oscillations be- 
low 14 T. 



II. EXPERIMENTAL DETAILS 

The single crystals of Cd-doped Bi2Se3 were grown by 
using elemental shots of Bi (99.9999%), Cd (99.99%) and 
Se (99.999%) as starting materials. To maximize the 
Cd content in the crystal, excess Cd and a mixture of 
Bi and Se with a ratio of Bi:Se — 32:68 were melted 
in a sealed evacuated quartz tube at 750 °C for 48 h 
with intermittent shaking to ensure a homogeneity, fol- 
lowed by cooling slowly to 550 °C and then annealing at 
the same temperature for one week. The resulting crys- 
tals are easily cleaved along the basal plane, revealing a 
silvery mirror-like surface. The X-ray diffraction mea- 
surements confirmed the crystal to be single phase with 
the proper Bi2Se3 structure. The actual Cd content was 
determined by the inductively-coupled plasma atomic- 
emission spectroscopy (ICP-AES) to be 0.0020(2). The 
as-grown crystals were examined by X-ray Laue analy- 



sis and cut into single-domain, thin bar-shaped samples 
with typical dimensions of 3x1x0.2 mm'^. 

For each annealing experiment, samples weighing 
about 4.2 mg were sealed in evacuated quartz tubes and 
annealed at a given temperature for one week, followed 
by quenching into cold water. All the samples used in this 
work were taken from the same part of the same batch, 
and the variation of the Cd content was confirmed to be 
negligible by the ICP-AES analysis. To avoid possible 
surface contamination, the surface layer of the annealed 
crystals were removed using adhesive tapes before trans- 
port measurements. 

It worth mentioning that in our annealing experi- 
ments, we took precautions to minimize the uncertainty 
in the annealing temperature. For each annealing run, 
we placed the quartz tube at the same position of the 
same furnace so that the temperature gradients in the 
furnace as well as the thermocouple calibration errors 
do not affect the annealing result. Also, the environment 
temperature was kept constant during this experiment to 
minimize the temperature fluctuations between different 
annealing runs. As a result, the annealing temperature 
Tannoai was vcry reproduciblc and its variation between 
different runs with nominally the same Tannoai was within 
±1 °C. 

The in-plane resistivity p^x and the Hall coefficient 
i?H were measured in a Quantum Design Physical Prop- 
erties Measurement System (PPMS-9) down to 1.8 K, 
for which the electrical contacts were prepared by using 
room-temperature-cured silver paste. In addition, one of 
the high-resistivity samples was brought to a 14-T mag- 
net for detailed SdH-oscillation measurements using an 
ac six-probe method, in which four lock-in amplifiers were 
employed to measure both the primary and the second- 
harmonic signals in the longitudinal and transverse chan- 
nels at a frequency of 19 Hz. The SdH-oscillation data 
were taken by sweeping the magnetic fields between ±14 
T with the rate of 0.3 T/min, during which the temper- 
ature was stabilized to within ±5 niK. 



III. RESULTS AND DISCUSSIONS 

A. p-type Bi2-a;Cdj;Se3 

The temperature dependence of the in-plane resistiv- 
ity Pxx of a Bii.g98Cdo.oo2Se3 crystal grown in the Se- 
rich condition is shown in Fig. 1(a), together with the 
data for a pristine Bi2Se3 sample grown with the same 
Bi/Se ratio. The pristine Bi2Se3 crystal shows an es- 
sentially metallic behavior with a weak resistivity up- 
turn below '^30 K, which is typical for low-carricr-density 
Bi2Se3;^^"^^ indeed, the Hall coefficient i?H in this sam- 
ple at 1.8 K corresponds to the bulk electron density 
no~7xl0^^ cm~'^, which is very small for Bi2Se3. On the 
other hand, in Bii.gg8Cdo.oo2Se3 the resistivity upturn is 
absent and the pxx value is lower, suggesting a higher 
carrier density. In fact, as shown in Fig. 1(b), i?H in the 
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FIG. 1: (Color online) (a) Temperature dependences of p^^^ 
of an as-grown Bii.998Cdo.oo2Se3 crystal (solid line) and a 
pristine Bi2Se3 crystal (dash-dotted line) grown in the same 
Se-rich condition. The pristine sample is n-type with n^ ~ 
7x10^^^ cm~^, while the Cd-doped sample is p-type with jih ~ 
9x10^* cm"''. Note that the low-temperature resistivity up- 
turn in the pristine sample is absent in the Cd-doped sample, 
(b) Temperature dependences of _Rh (left axis) and ^h for the 
Cd-doped sample; inset shows the magnetic field dependence 
of py:c in this sample at 1.8 K. 



Bii.9g8Cdo.oo2Se3 crystal is positive and its value at 1.8 
K corresponds to the hole density nh~9xl0^* cm~'^, im- 
plying that the Cd doping has created ~9.7xl0^^ cm~^ 
of holes (^^0.8 hole per Cd atom) that outnumbers the 
electrons coming from Se vacancies J^ The reason for the 
smaller number of doped holes compared to the Cd con- 
centration is most likely that a small portion of Cd atoms 
occupy the interstitial site and act as donors. In Fig. 
1(b), i?H is only weakly dependent on temperature and 
the Hall resistivity pyx is perfectly linear in B (as shown 
in the inset), refiecting the metallic nature of the as- 
grown Bii.gg8Cdo.oo2Se3 sample which is due to a single 
type of carriers (i.e., the bulk holes). The Hall mobility 
/iH (= Rn/ Pxx), also shown in Fig. 1(b), increases with 
decreasing temperature, reaching ~1600 cm^/Vs at 1.8 
K. 



B. p-to-n-type conversion by post annealing 

Annealing the as-grown BiigggCdo. 002863 crystals in 
evacuated quartz tubes has a drastic effect on its trans- 
port properties. Figure 2(a) shows how the tempera- 
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FIG. 2: (Color online) (a) Temperature dependences of pxx 
for Bii,9g8Cdo.oo2Se3 crystals annealed at different tempera- 
tures in evacuated quartz tubes. The data for the as-grown 
crystal is also shown for comparison. The low-temperature 
resistivity values span three orders of magnitude, (b) Tem- 
perature dependences of Rh for the same series of samples. 
As the annealing temperature is increased, the crystal be- 
comes less p-type and is eventually converted to n-type. This 
conversion occurs with the annealing temperature of ^--^577 
°C. 



ture dependence of pxx changes upon annealing in a nar- 
row temperature window between 573 and 590 °C. One 
can see that Pxx{T) evolves nonmonotonically with the 
annealing temperature, Tanneai; namely, pxx initially in- 
creases with Tannoai Until Tanneai cxcccds 577 °C, after 
which Pxx decreases as Tanneai is further increased. No- 
tably, Pxx of the sample annealed at 577 °C show a high 
value of 0.5 17cm at 1.8 K, which is three orders of magni- 
tude larger than that of the as-grown sample. Figure 2(b) 
shows the typical temperature dependences of low-field 
i?H data (defined as i?H — Pyx/B for B ^ 0) for differ- 
ent Tannoai J wliich indicates that more and more electron 
carriers are introduced as Tanneai is increased, and the 
sign change from p-type to n-type occurs around Tannoai 
= 577 °C. 

Since the drastic change in the transport properties 
occurs in a very narrow temperature window (573 - 580 
°C), one may wonder about the reproducibility of the 
result. As a matter of fact, the observed change was 
quite reproducible, as demonstrated in Figs. 3 and 4. 
Figures 3(a)-3(d) show the Pxx[T) data for at least two 
samples annealed at the same temperature, where one 
can see that the behavior for each Tannoai is essentially 
reproducible. Figures 4(a)-4(d) show the corresponding 
i?H(T) data for the same sets of samples; here, except 
for the case of Tannoai — 577 °C [Fig. 4(c)], we observed 
reasonable reproducibility [Figs. 4(a), 4(b), and 4(d)]. 
For Tannoai — 577 °C, two of the three samples (A and 
B) showed a sign change in i?H from negative to positive 
upon lowering temperature from 300 K, whereas the i?H 
of sample C remained negative in the whole temperature 
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FIG. 3: (Color online) (a-d) Reproducibility of the pxx(T) 
data in samples annealed at the same temperature, demon- 
strated for four different values of Tanneai indicated in each 
panel. 



range. Actually, this variation in the behavior of i?H 
indicates that the Tannoai = 577 °C samples are at the 
verge of the p-to-n-type conversion. 

Figures 4(e)-4(h) show the Pyx{B) curves measured in 
sample "A" of each Tanneai- One can clearly see that the 
curve in Fig. 4(g) for Tannoai = 577 °C is nonlinear, in- 
dicating that there are at least two bands contributing 
to the transport. In topological-insulator samples with a 
large bulk resistivity, this kind of nonlinear pyx{B) curves 
are indications of the surface channels making noticeable 
contributions. ^^"^^ Therefore, the consistently high resis- 
tivity [Fig. 3(c)] together with the complex behaviors 
of the Hall signal are likely to be a signature of a high 
level of compensation achieved in the samples annealed 
at 577 °C; in other words, in those samples the acceptors 
and donors are nearly equal in number and their delicate 
balance can easily change the sign of R-a ■ It is to be em- 
phasized that our data demonstrate that this high level 
of compensation is reproducibly achieved with Tannoai = 
577 °C. In samples annealed at other temperatures, the 
Pyx{B) behavior is almost linear [Figs. 4(e), 4(f), and 
4(h)], suggesting that the contribution of the surface to 
the transport properties is minor. In passing, the col- 
lection of Pxx{T) and Rk{T) data shown in Fig. 2 for 
varying Tanneai are for sample "A" of each Tannoai shown 
in Figs. 3 and 4. 



C. Defect chemistry 

The above observation that a drastic change in the 
transport properties of Bii.998Cdo.oo2Se3 occurs in a 
very narrow temperature window might seem surpris- 
ing. However, this behavior can be readily understood 
by examining the defect chemistry associated with the 
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FIG. 4: (Color online) (a-d) Reproducibility of the Ru{T) data in samples annealed at the same temperature, demonstrated 
for four different values of T^^nncai indicated in each panel, (e-f ) Magnetic-field dependences of pyx measured in sample "A" of 
each Tanncai- The solid line in (g) is the two-band-model fitting to the data. 



annealing. In the present system, there are mainly two 
different types of charged defects, the aliovalent substi- 



tutional defect Cdgj^ and the vacancy defect V** f 



■ the 



former acts as an acceptor and the latter as a donor. 
Therefore, the effective charge-carrier density is deter- 
mined by their competition and can be expressed as rioff 
— [Cdgi] — 2[V^"], where positive (negative) ricff' denotes 
the hole (electron) density. In an as-grown sample, the 
Cdgi defects are dominant and UeS is positive; accord- 
ingly, the chemical potential lies in the valence band. 
When annealed in evacuated quartz tubes, a portion of 
selenium goes into the gas phase Se2 in equilibrium with 
the solid phased resulting in the formation of more V** 
defects while leaving Cdg^^ unaffected (because the Ta„ncai 
employed in the present study is much lower than the 
melting temperature of 710 °C). The equilibrated vapor 
pressure of Se2 increases with increasing Ti^nncaii^'^ creat- 
ing more V^" and eventually changing rics from positive 
to negative. 

To be more quantitative, one may assume that the in- 
crease in the Se-vacancy concentration upon annealing, 
A[l^**], is directly reflected in the increase in the number 
of Se2 molecules in the quartz tube, which determines the 
Se2 vapor pressure Psca! in constant volume, one expects 
a linear relation between A[Vse*] and Pso2 if the Sep va- 
por behaves as an ideal gas. According to Ref. [2j, the 
equilibrated Se2 vapor pressure Pso2 of Bi2Se3 is related 
to the absolute temperature T via 



logPseJatm]=A-B/r[K], 



(1) 



where A = 7.81±0.50 and B = 10870±640 for the tem- 
perature range of 527 to 627 °C. From this T dependence 
of Psc2: one can infer that A [14"] is very sensitive to the 



change in Tanncai: For example, changing Tanneai by just 1 
°C near 577 °C results in a variation of ~3.5% in A[V^s"]; 
therefore, the expected change in A[V^"] upon changing 
Tanncai from 575 to 580 °C is as much as ^18%. It is thus 
expected that a sign change in noff occurs abruptly in 
the vicinity of Tanncai = 577 °C where UcS ~ (namely, 
[Cdg^] ~ [V"]), and results in a drastic change in the 
transport properties as we observed. 



D. Surface quantum oscillations 

From the above results, it is clear that the highest level 
of compensation is achieved in samples annealed at 577 
°C. We therefore measured the sample C of Tanncai — 577 
°C in a 14-T magnet using a rotation sample holder to 
investigate its SdH oscillations in detail. Before the high- 
field measurements, to protect the surface state from ag- 
ing, the top surface of the sample was covered with AI2O3 
in the following way: first, the crystal was cleaved on 
both surfaces with adhesive tapes to reveal fresh sur- 
faces, mounted on a sample holder with GE varnish, and 
transferred into the sputtering chamber; second, the top 
surface was cleaned by bias-sputtering with Ar ions for 13 
minutes and then, without breaking the vacuum, a 540- 
nm-thick AI2O3 film was deposited by the rf magnetron 
sputtering. After this process, gold wires were bounded 
to the side faces by spot welding. Probably because of the 
sample heating during the spot welding, the pxx value of 
this sample became even larger than that shown in Fig. 
3(c). Also, the sign of Ru at low temperature changed 
to negative after the process, showing the Hall response 
similar to that of the sample A. The i?H (T) behavior of 
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FIG. 5: (Color online) (a) Rh{T) data of the sample C of 
Tknneai = 577 °C after the AbOa-coverage process, (b) Pyx{B) 
data of the same sample measured at 1.4 K; the solid line 
is the result of the two-band-model fitting, (c) Magnetic- 
field dependences of the primary signal {pxx) and the second- 
harmonic signal [Apxx) measured at 1.4 K in magnetic fields 
along the C3 axis, (d) Plots of l^pxx{B) and d?pxx/dB^ vs 
the inverse magnetic field 1/ B. Clear SdH oscillations can be 
seen in both data, and good agreements in the positions of 
peaks and dips between the two curves are evident, (e) The 
FT spectrum of /S.pxx showing three prominent frequencies. 



this sample C after the Al2 03-coverage process is shown 
in Fig. 5(a), and its pyx{B) curve at 1.4 K is shown in 
Fig. 5(b). 

In our measurements using an ac lock-in technique, 
we simultaneously recorded the primary and the second- 
harmonic signals during the magnetic-field sweeps. Fig- 
ure 5(c) shows the magnetic-field {B) dependences of the 
primary signal {pxx) and the second- harmonic signal (de- 
noted /S.pxx) measured at 1.4 K in magnetic fields along 
the C3 axis. One can see that the second-harmonic signal 
{l^Pxx) shows pronounced oscillations, while in the pri- 
mary signal {pxx) the oscillations are hardly visible. To 
understand the nature of the oscillations, we show in Fig. 
5(d) the plots oi I^pxx and d^pxx/dB^ (second derivative 
of the primary signal) vs the inverse magnetic field 1/B; 
a comparison between the two curves indicates that they 
present essentially the same peak/dip positions. While 
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FIG. 6: (Color online) (a) The Apxx data for varying 
magnetic-field directions, plotted as a function of 1/(B cos6); 
dashed lines mark the positions of the peaks, (b) The three 
prominent frequencies in the FT spectra of the SdH oscilla- 
tions plotted as a function of 9. All the frequencies vary as 
1/ cos 6, as indicated by the solid lines. 



the waveforms are quite complicated, the Fourier trans- 
form (FT) spectrum of Apxx{B^^) shown in Fig. 5(e) 
presents three well-defined peaks at Fi = 11 T, F2 = 27 
T, and F3 = 33 T, indicating that the observed oscilla- 
tions are SdH oscillations with multiple frequencies. 

Note that the second-harmonic in ac measurements is 
a distortion of the input sine wave, and its occurrence 
is an indication of a nonlinear response. In the present 
case, the SdH oscillations are apparently giving rise to a 
peculiar non-ohmicity. This makes the second-harmonic 
signal to be useful for observing the SdH oscillations with 
a high sensitivity, although the detailed mechanism is not 
clear at the moment. 

Figure 6(a) shows how the SdH oscillations observed 
in Apxx change when the magnetic field is rotated, by 
plotting Apxx versus 1/{B cos9) where 6 is the angle be- 
tween B and the C3 axis. One can see that the oscillatory 
features are essentially dependent on the perpendicular 
component of the magnetic field. Also, as shown in Fig. 
6(b), the angular dependences of all three frequencies in 
the FT spectra are consistent with 1/ cos 6*. These results 
strongly suggest that the present SdH oscillations signify 
2D Fermi surface(s). We note that the SdH oscillations in 
this sample disappeared after keeping the sample in am- 
bient atmosphere for a week, which suggests that the SdH 
oscillations were coming from the surface. (This observa- 
tion also suggests that the AI2O3 coverage, while useful 
for slowing the aging of the surface of Bi2Se3, does not 
provide a perfect protection). Furthermore, one can esti- 
mate the bulk mobility of this sample to be ~40 cm^/Vs 
from the values of pxx and Rn , and such a mobility is too 
low to give rise to SdH oscillations of the bulk carriers 
below 14 T. All told, one can reasonably conclude that 
the observed SdH oscillations are of the surface origin. 



Prom the SdH-oscillation data, the Fermi wave vec- 
tor kp can be calculated via the Onsager relation F — 
ihc/2Tre)TTkl, yielding kp = 0.018, 0.029, and 0.032 A"! 
for Fi, F2, and F3, respectively. These are of the same 
order as the value kp — 0.031 A~^ reported for the topo- 
logical surface state in n-type Bi2-2;Sba:Se3.— However, 
because of the multi-component nature of the oscillations 
that leads to complicated waveforms, it is difficult to re- 
liably extract the cyclotron mass nic nor the Dingle tem- 
perature for each component using the Lifshitz-Kosevich 
theory.^'* This makes it impossible to identify the origins 
of the three oscillation frequencies, but possible reasons 
for the multiple components in the present case include: 
(i) harmonics of a fundamental frequency are observed, 
(ii) chemical potentials of the top and bottom surfaces 
are not identical and give two frequencies associated with 
the topological surface states, (iii) a trivial 2D electron 
gas^^ created by the band bending at the surface presents 
additional SdH oscillations. Note that, in the case of the 
present sample used for the detailed SdH measurements, 
the top surface was covered by AI2O3 and the bottom 
surface was covered by the GE varnish, so the conditions 
of the two surfaces were very different. To resolve the 
origins of those multiple frequencies, an experiment in- 
volving the gate control of the surface chemical potential 
to trace the energy dispersion of each branch (as was 
done^ for exfoliated Bi2Se3) would be desirable. 



E. Nonlinear pyx{B) behavior 

Although we could not extract the surface mobility 
from the SdH oscillations in the present case, we can 
still estimate the relevant parameters of the surface and 
the bulk transport channels by analyzing the nonlinear 
Pyx{B) behavior [Fig. 5(b)] with the simple two-band 
model described in Ref. 18. The solid line in Fig. 5(b) is 
the result of the two-band-model fitting, from which we 
obtained the bulk electron density rib — 7x10^^ cm~'^, 
the bulk mobility /i;, = 17 cm^/Vs, the surface electron 
density Us — 1.8x10^^ cm~^, and the surface mobility 
/is = 1.2x10^^ cm^/Vs. The observed frequencies of the 
SdH oscillations, 11, 27, and 33 T correspond to the sur- 



face carrier densities of 2.5x10^^, 6.7x10^^, and 8.1x10^^ 
cm~^ in spin-filtered surface states, respectively, and it 
is interesting that the sum of these numbers, 1.7x10^^ 
cm~^, appears to be consistent with the Ug value ob- 
tained from the two-band analysis. Also, it is assuring 
that the surface mobility fig obtained from the two-band 
analysis, ■~1200 cm^/Vs, is reasonably large and is con- 
sistent with our observation of the SdH oscillations in 
moderate magnetic fields. 

In passing, the pyx{B) data of the sample A shown in 
Fig. 4(g) can also be fitted with the two-band model. 
The result of the fitting, shown by the solid line in Fig. 
4(g), yields the bulk electron density rif, = 1.3x10^* 
cm"'^, the bulk mobility pb = 15 cm^/Vs, the surface 
electron density Us = 2.0x10^"'^ cm~^, and the surface 
mobility /i^ = 1.0x10^ cm^/Vs for this sample. 



IV. CONCLUSION 

In conclusion, we demonstrate that with tactful de- 
fect engineering one can optimize the transport proper- 
ties of the topological insulator Bi2Se3 to show a high 
bulk resistivity and clear quantum oscillations. Specifi- 
cally, by employing a Se-rich crystal-growth condition we 
achieved the p-type state in Bi2Se3 by Cd-doping for the 
first time; we then employed careful post annealing to 
tune the Se vacancies and achieved a high level of com- 
pensation, where the acceptors and donors nearly cancel 
each other and the sample presents a high p^x value ex- 
ceeding 0.5 J7cm at 1.8 K and shows 2D SdH oscillations 
consisting of multiple components below 14 T. 
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